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ABSTRACT 

A bending fracture test specimen design is presented for 
thin elements used in optical devices for space flight appli- 
cations. The specimen design is insensitive to load position, 
avoids end effect complications, and can be used to measure 
strength of membranes less than 2 /xm thick. The theoreti- 
cal equations predicting stress at failure are presented, and 
a detailed finite element model is developed to validate the 
equations for this application. An experimental procedure us- 
ing a focused ion beam machine is outlined, and results from 
preliminary tests of 1.9 /xm thick single crystal silicon are pre- 
sented. These tests are placed in the context of a methodol- 
ogy for the design and evaluation of mission critical devices 
comprised of large arrays of cells. 


NOMENCLATURE 



Figure 1: Micro-shutter element suspended from torsion 
bars, shown rotated almost 180°. 
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1 MOTIVATION 

NASA/Goddard Space Flight Center develops Micro- 
electromechanical System (MEMS) optical components for 
cryogenic space applications. Many of the devices are etched 
from 2.0 /xm thick single crystal silicon or 0.5 /xm thick silicon 


nitride membranes, involve extremely large deformations, and 
sustain huge permanent or periodic stresses during cooling 
and operation. 

Examples of these applications are shown in Figs. 1 and 2. 
The first case is a large micro-shutter array patterned from 
0.5 /xm thick silicon nitride with a one side Aluminum layer of 
about 20 nm [1 ' 2] . The single shutter shown in Fig. 1 is com- 
prised of a 100 /xm X 100 \im blade suspended from a 2.0 /xm 
wide torsion bar fixed at the ends. During operation the shut- 
ter is cooled to 30 K and repeatedly rotated 90°. The shutter 
is shown rotated almost 180° in the above figure, illustrating 
the flexibility of this material at these dimensions. Detailed fi- 
nite element models based on bulk properties of silicon nitride 
predict a peak stress of 4 GPa t yet the bulk strength of hot 
pressed silicon nitride is on the order of 800 MPa. 

The second example shown in Fig. 2 is an array of “pop-up” 



Figure 2: Pop-up Detectors shown in the permanent 
"folded'* configuration. 



Figure 3: Dual width bending test specimen for brittle 
membranes less than 2 fim thick. 


detector elements^. The suspended pixel and support leg 
structures are initially etched from a flat single crystalline sil- 
icon membrane of 2 fim nominal thickness. Each element 
contains an embedded bolometer. The central pixel area is 
coated with a 10 jam thick layer of absorbing material such 
as bismuth. In the last fabrication step, the support legs are 
folded 90°, resulting in the free standing structures shown in 
Fig. 2. During operation the detectors are cooled to less than 
IK. A detailed finite element model predicts a peak stress of 
2 GPa at the stress concentration due to the folding process, 
yet the bulk strength of single crystal silicon is 200-300 MPa. 

Each of these device designs exploit the volume effect on 
strength of brittle materials. Initial design of the devices has 
been completed using trial and error, and qualitative sim- 
ulations based on conservative estimates of stiffness and 
strength. Quantitative simulations based on accurate stiffness 
and statistical strength data are required to estimate the prob- 
ability of device failure and to help guarantee device reliability 
during space flight. The mechanical material properties, in 
particular strength, vary with size and processing. As a result, 
inexpensive and convenient methods to measure properties 
of a large numbers of samples are needed. 

In this paper we present the details of a dual width bending 
specimen used to measure the strength of thin membranes. 
Subsequently we discuss the application of this approach for 
the design and analysis of optical components etched from 
thin membranes. 


2 DUAL WIDTH FRACTURE TEST 

Simple cantilever beam test specimens suffer from the follow- 
ing limitations for thicknesses less than 2 jim\ 


1 . Applied force is too small to measure directly. 

2. Results strongly depend on force position. 

3. Fixed end failure is complicated by end effects. 

4. Very thin silicon and silicon nitride beams are too flexible 
and strong to break in bending at reasonable angles. 


The dual width bending specimen of thickness t shown in 
Fig. 3 is designed to overcome these limitations. At the fixed 
end, the beam has a width w Wt which necks down to narrow 
width w n for a length L n . The remaining portion of the beam 
returns to width w w . For a vertically applied load F at point 
d the peak stress will occur along the section at point 6 for 
w w /wn > L a /L b . Failure away from the fixed end avoids the 
complicated end effects mentioned in the third item above. In 
this configuration, the rotations at points b and c, 9 b and 0 C , re- 
spectively, are measured just prior to failure using image pro- 
cessing. These rotations along with the elastic properties are 
used to compute the peak stress at point b, without knowing 
force F. 

We also select w w » w n such that the wide portion of the 
beam shown in Fig. 3 is rigid relative to the neck region. This 
drives most deformation into the the neck region, resulting in 
higher stresses at lower deflections, alleviating the last item 
above. This configuration also reduces the sensitivity to force 
location as shown in the following section. 


2.1 Theoretical Equations 

We may compute the relative rotation 9 C - 9 b of the neck region 
from beam theory, and use this to compute the peak stress at 
section b. For a vertical force F applied through the rigid wide 
beam at a distance L b from point 6, the relative rotation is 



Undeformed 



Figure 4: Detailed finite element model one symmetric 
half of the dual width fracture test specimen. 
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Figure 5: Deformations predicted by the finite element 
model. 


where L n is the length of the neck region, E is the Young’s 
modulus, and I = w n t 3 /12 is the cross section area moment 
of inertia of the neck region. For Lb ~ aL n the above equation 
can be rewritten as 


(2) 


and pick an arbitrary failure stress 07 = 4 GPa , which is dou- 
ble the average value of 2 GPa presented in Ref. 4. We align 
the length of the specimen with the <1 10> direction and as- 
sume a modulus of E = 176 GPa . For a thickness t = 1.9 fim, 
inclination 0 C - 0 b - 10°, and scale factor a = 10, equation (4) 
yields L n = 7.5 pm. Therefore, a beam neck less than this 
length should result in failure at relatively small rotations. For 
convenience we chose L n — 6.55 fim for the remaining stud- 
ies in this paper. 


The resulting maximum moment M b about point b is 


M b = FL b 


EljOc - Ob) 
Ln (1 -&)' 


The above equations are derived based on linear beam the- 
ory. A detailed plate finite element model is presented in the 
next section to validate the above equations for the dual width 
specimen configuration under large deflections, using the pa- 
rameters specified above. Subsequently, test results are pre- 
sented and used to compute the strength. 


The peak stress occurs at section b and is = Mbt/2I. 
Substituting equation (3) results in 


eq _ Et(0 c Ob) 
XX ~ 2L n ( 1 -&) 


( 4 ) 


The above equation is insensitive to force position for large 
a = L b /Ln . Even for smaller a > 2, a small error in a has lit- 
tle effect on the results. Consequently, only approximate mea- 
surement of the distance L b is required to obtain reasonably 
accurate measures of stress. Small variations in force position 
during deformation will also have little effect on the results. 

Equation (4) can be used to choose L n such that the speci- 
men will fail for small rotations. For the purpose of designing 
the specimen we consider a single crystal silicon membrane 


2.2 Finite Element Model 

A detailed finite element model of the dual width specimen 
with parameters given above, widths w n = 6.25 fim and 
w w = 30 fim, and inside corner radius of 1 /im is shown 
in Fig. 4. Only one symmetric half of the specimen geome- 
try is modeled to reduce computation time and improve non- 
linear solution convergence. The finite element model is con- 
structed of 166 four-node plate elements, and solutions were 
obtained using the non-linear geometry solution sequence in 
UAI/NASTRAN [5] . 

For the orthotropic material properties, we use moduli E = 
176 GPa and G - 80 GPa corresponding to the the <1 10> 
direction, moduli E = 130 GPa and G = 50 GPa correspond- 
ing to the <100> direction, and Poisson’s ratios of v = 0.279 
and v = 0.064, respectively. In this case the beam is located 
in the {100} plane aligned with the <11 0> direction. All 6 de- 
grees of freedom are fixed along the left edge and a normal 
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Figure 6: Bending test specimen used to measure strength. 


pressure is applied to one element as shown in Fig. 4. This 
load is treated as a follower force such that the load remains 
normal to the element surface during deformation. 

Deformation results predicted by the finite element model at 
a relative rotation of 0 C - 0 b = 11° are shown at true scale in 
Fig. 5. The above figure illustrates that the end deflection is 
relatively large, and that slight rotation does occur at point b as 
expected. Stress contour fringes of the finite element stress 
predictions in the x direction are shown in the neck region 
in Fig. 6. As expected, the peak stress occurs in the neck 
region. 

The maximum principal stress and x component of stress 
ah are listed in Table 1 for increasing relative rotation 6 C - 0 b . 
Corresponding values of <r% q x computed from equation (4) are 
also listed and compared to the finite element predictions. The 
theoretical and finite element results are in excellent agree- 
ment for <t xz . Brittle failure will most likely be driven by the 
maximum principal stress which is about 10% higher than the 
x component. As a result the theoretical equation will under- 
predict strength. This is acceptable and even desirable be- 
cause it is conservative. 

There is a concentration in stress due to the inside corner. The 
finite element method underpredicts the peak stresses due to 
averaging over each element. Equation (4) does not account 
for the stress concentration and as a result yeilds conserva- 
tive estimates of strength. Hence there is good agreement 
through cancellation of errors. As a result the strength pre- 
dictions presented herein are well suited for comparison with 
finite element results for actual structures such as the micro- 
shutters and pop-up detectors. 

The above analysis validates equation (4) for predicting the 
strength using the dual width fracture test specimen. In the 
next section we use this equation to estimate the strength of 
single crystal silicon using the fracture test. 


TABLE 1: Comparison of results for w n — 6.25 tim, 
L n = 6.55 jJim, Wu, = 30 fim, L a = 73.5 /im. 
Lb = 58.5 ixm , a = 8.9. 


2.3 Experimental Measurements 

All experiments were performed in an FEI 620 focused ion 
beam milling machine. The FEI 620 is a dual beam machine, 
with an ion and electron column, which allows ion milling and 
in situ scanning electron microscopy. Ion milling can be per- 
formed with ion beam spot sizes of 20 nm to 1 //m and ma- 
chining rates in silicon up to 10 n m 3 /s . The machine is also 
equipped with a micro-manipulator needle and the ability to 
deposit platinum by ion induced metal organic chemical vapor 
deposition (MOCVD). 

The above platform was used to machine and test the dual 
width geometry from single crystal silicon membranes with 
the parameters specified in the previous section. The micro- 
manipulator is used to gradually apply the vertical load F to 
failure. A series of scanning electron microscope (SEM) im- 
ages are recorded, and the last image prior to fracture is an- 
alyzed to extract the rotation angles 6 C and 6 b . These angles 
are used in conjunction with equation (4) to compute the peak 
stress. 

Prior to bending, a light line is etched through the center of the 
neck region using the focused ion beam (shown as a dashed 
line in Fig. 3). A deformed cantilever specimen is shown in 
Fig. 7. The image just prior to failure is analyzed by drawing 
lines tangent to the etch line at points b and c as shown in 
Fig. 7. Horizontal and vertical reference lines are added. The 
base and height dimensions of the resulting triangles are mea- 
sured and transformed to account for the image angle relative 
to the specimen plane. These dimensions are used to calcu- 
late the rotation angles. 

Using the above procedure, the relative rotations were mea- 
sured just prior to failure for 9 specimens, and the correspond- 
ing peak stress was calculated. The results are listed in Ta- 






Figure 7: Scanning electron micrograph of a bending test 
specimen used to measure strength. 


ble 2. The average stress at failure was found to be 14.2 GPa 
with a standard deviation of 1.37 GPa. 

Wilson et. al. 3 [4] presented strengths in the range of 1-2 GPa 
for micro-cantilever silicon beams with thickness of 30 fim. 
However, in the study presented herein, the membrane thick- 
ness is 1.9 /im, and the membranes are coated with metallic 
substances. These aspects will lead to strengths that are dif- 
ferent from values presented elsewhere. The membranes are 
also so thin that simple cantilever beam tests tend to bend be- 
yond 90° before failure, making it difficult to ascertain strength. 

The above failure stress values are on the same order as 
those reported by Johansson et. alJ 6] , but much greater than 
the values presented by Wilson et. al. [4} . This is consistent 
with the fracture theory that smaller volume results in fewer 
flaws and higher strengths. Here we have thicknesses less 
than one tenth of the thicknesses used by Wilson et. al. [4] . 
The coating material may also provide a small increase of the 
strength t6] . 

Here we note that the results presented in Table 2 are highly 
approximate. In practice we found it difficult to accurately mea- 
sure the rotation angles directly. As a result the rotations in Ta- 
ble 2 are only within +/- 5 degrees, which represents an error 
of +/- 1 5 to 20% for the angles listed in the table. The stresses 
computed using equation (4) are directly proportional to these 
angles, which means the listed stresses are within +/- 15 to 
20%. Ultimately, this method requires improved methods of 
determining the rotation angles to yield reliable results. 



TABLE 2: Fracture stresses computed from experimental 
measurements f or t = 1.9 /um. The average is 14.2 GPa 
and the standard deviation is 1.37 GPa. The angles were 
measured with an accuracy of +/- 5 degrees. 


data in material libraries, handbooks, or archived literature. 
We typically use such cataloged data for the purposes of de- 
sign of large scale brittle structures, perhaps complimented 
with some critical testing prior to fabrication. Once the struc- 
ture is fabricated we inspect and test a prototype or produc- 
tion unit under environments that envelope anticipated service 
loads. 

There are key drawbacks to this procedure for the develop- 
ment of MEMS devices that endure substantial stresses. Most 
importantly, strength is strongly a function of processing. In 
large scale structures processing effects are minimized by 
eliminating edges from the design and applying post fabrica- 
tion surface treatments to reduce the population of surface 
flaws. MEMS edges left by etching can not be removed, and 
there are limited post etching treatments available. MEMS de- 
vices are also extremely sensitive to process variations across 
a wafer processed in batch. For example, the micro-shutter 
array described earlier will ultimately be comprised of at least 
1024 by 1024 elements. At 100 \im each this adds up to about 
10 cm. As a result, for applications where reliable performance 
is critical the above traditional procedure is not adequate. 

Development of such devices requires simple and convenient 
tests, even at the expense of some accuracy, as long as the 
results are conservative. Initially a large number of test sam- 
ples may be fabricated and tested with the specific process 
anticipated for a given device. This data may be used to com- 
plete the design using Weibull statistical analysis as described 
in Ref. 7. Ultimately, a few witness samples could then be 
included at various locations in the actual device fabrication. 
Tests from these witness samples can then be compared with 
the larger data set to validate the resulting process and ulti- 
mately the actual strength of the fabricated device. 


3 FRACTURE PROBABILITY METHODOLOGY 

Traditionally researchers use the most accurate test methods 
available to measure and catalogue statistical brittle strength 



4 CONCLUSION 

In this work we presented a dual width fracture test specimen 
that helps alleviate several difficulties associated with simple 
cantilever tests. In particular, toad position sensitivity is re- 
duced and complications due to end effects are eliminated. 
This specimen geometry also allows the strength to be ap- 
proximated without knowing loads which are difficult to re- 
solve. This specimen may be used for membranes less than 
2 jim thick. Bending test methods are attractive because they 
provide relatively convenient witness samples. 

While this specimen does reduce the force position sensitiv- 
ity needed to accurately compute stress, the test does suffer 
from some of the same shortcomings of any bending test at 
these scales. In particular, the results are strongly dependent 
on accurate metrology, and are based on a presumed elastic 
modulus. There are also some concerns regarding the merits 
of bending tests compared to uniaxial tension tests [8 l Failure 
in bending is due to highly localized stresses, and does not 
account for the overall distribution of flaws in a volume, which 
may yield overly optimistic results. In this case, however, the 
bending test closely mimics actual load configurations in our 
devices. 

The results presented here are promising for future investiga- 
tions. In particular, each of our devices operate at cryogenic 
temperatures, and it will be important to understand the tem- 
perature effects on strength. These tests are also amenable 
to fatigue testing by placing an electrode near the specimen 
and cycling the applied voltage. These results are inherently 
less accurate than uniaxial tension tests which measure both 
force and deflection [9] . As a result, it will be important to rig- 
orously validate these results with more accurate uniaxial ten- 
sion tests. 
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